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ABSTRACT 

Context. AB Doradus is the main system of the AB Doradus moving group. It is a quadruple system formed by two widely separated 
binaries of pre-main-sequence (PMS) stars: AB Dor A/C and ABDorBa/Bb. The pair AB Dor A/C has been extensively studied and 
its dynamical masses have been determined with high precision, thus making of AB Dor C a benchmark for calibrating PMS stellar 
models. If the orbit and dynamical masses of the pair AB Dor Ba/Bb can be determined, they could not only play a similar role to that 
of AB Dor C in calibrating PMS models, but would also help to better understand the dynamics of the whole AB Doradus system. 
Aims. We aim to determine the individual masses of the pair AB Dor Ba/Bb using VLBI observations and archive infrared data, as 
part of a larger program directed to monitor binary systems in the AB Doradus moving group. 

Methods. We observed the system AB DorB between 2007 and 2013 with the Australian Long Baseline Array (LBA), at a frequency 
of 8.4 GHz in phase-reference mode. 

Results. We detected, for the first time, compact radio emission from both stars in the binary, AB DorBa and AB DorBb. This result 
allowed us to determine the orbital parameters of both the relative and absolute orbits and, consequently, their individual dynamical 
masses: 0.28±0.05 M 0 and 0.25±0.05 M 0 , respectively. 

Conclusions. Comparisons of the dynamical masses with the prediction of PMS evolutionary models show that the models underpre¬ 
dict the dynamical masses of the binary components Ba and Bb by 10-30% and 10-40%, respectively, although they all still agree 
at the 2 -a level. Some of stellar models considered favour an age between 50 and 100 Myr for this system, meanwhile others predict 
older ages. We also discuss the evolutionary status of ABDorBa/Bb in terms of an earlier double-double star scenario that might 
explain the strong radio emission detected in both components. 

Key words, stars: individual: AB DorB - binary stars - radio emission 


1. Introduction 

Stellar evolution models are used to predict fundamental pa¬ 
rameters of the stars, such as their mass. The estimates from 
mass-luminosity theoretical relationships are not in concordance 
for the particular case of pre-main-sequence (PMS) stars with 
masses < 1.2 M e . Either more accurate observations or re¬ 
vised models are needed. As reported in Azulay et al. (2014), 
a VLA/VLBI program to detect binary stars with substantial 
emission at radio wavelengths is underway. This program is 
focused on stars members of the AB Doradus moving group, 
AB Dor-MG (Zuckerman et al. 2004), which includes the al¬ 
ready studied pairs AB Dor A/C (Guirado et al. 2006, G06) and 
HD 160934 A/c (Azulay et al. 2014). 

Located at a distance of ~15pc (G06), the stellar system 
AB Doradus has a pair of binaries separated by ~9", AB Dor A 
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and ABDorB. The K1 star ABDorA, the main star of the 
system, is an active star with detectable emission at all wave¬ 
lengths. In particular, it has strong radio emission (Guirado et al. 
1997) generated by a dynamo effect caused by its rapid rotation 
(0.514 days). This star has a low-mass companion, ABDorC 
(0.090Mo), whose study is important to calibrate the stellar evo¬ 
lution models of young low-mass stars (Close et al. 2005). 

The other binary of the AB Doradus system is ABDorB 
(=Rossiter 137 B), consisting of two components, AB DorBa 
and AB DorBb, with spectral types M5 and M5-6 (Close et al. 
2007), respectively. These components are separated by an angu¬ 
lar distance of ~0.06" (G06; Janson et al. 2007). The combined 
system shows a high rotation rate with a period of 0.38 days 
(Lim 1993; Wolter et al. 2014) and displays strong radio emis¬ 
sion detected both by the Australian Telescope Compact Array 
(ATCA; Lim 1993; Wolter et al. 2014) and the Australian VLBI 
Network (G06; this paper). Close et al. 2007 and Wolter et al. 
(2014) report information on the relative orbit from NIR VLT 
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Table 1. Journal of observations of AB DorB 


Date (Epoch) 

Array configuration 11 

UT Range 

Beam size 
[mas] 

PA. 

[°] 

11 nov 2007 (2007.863) 

At, Cd, Ho. Mp, Pa, Hh 

10:00-22:00 

3.04x1.18 

-2.4 

25 oct 2010 (2010.816) 

At, Cd, Ho. Mp, Pa 

11:00-23:00 

2.96x2.76 

74.3 

16 aug 2013 (2013.625) 

At, Cd, Ho, Mp, Pa, Hh, Ti, Ww 

15:00-03:00 

2.61x1.05 

0.8 


Notes. a : At: Australia Telescope Compact Array, Cd: Ceduna, Ho: Hobart, Mp: Mopra, Pa: Parkes, Ti: DSS43 - NASAs Deep Space Network 
Tidbinbilla, Ww: Warkworth, Hh: Hartebeesthoek 


observations at H and K bands. The latter authors estimated a 
value of the sum of the masses of the components Ba and Bb 
of 0.69^ ^ M G , somewhat larger than the model-dependent es¬ 
timates from Janson et al. (2007) (0.13-0.20 and 0.11-0.18 M G 
for Ba and Bb, respectively). 

An important point of debate is the age of the AB Doradus 
system, which is not determined with sufficient accuracy and 
presents different estimates in different publications: from 40- 
50Myr (Zuckerman et al. 2004; Lpez-Santiago et al. 2006; 
Guirado et al. 2011) to 50-100Myr (Nielsen et al. 2005; Janson 
et al. 2007; Boccaletti et al. 2008) and 100-140Myr (Luhman et 
al. 2005; Ortega et al. 2007; Barenfeld et al. 2013). The deter¬ 
mination of the age of the AB Doradus system is fundamental to 
calibrate the evolution models of PMS stars. 

In this paper, we report the results of three epochs of VLBI 
observations of ABDorB, leading to the discovery of radio 
emission from both components, the determination of the dy¬ 
namical masses for the individual components, and their com¬ 
parisons with theoretical models. We also discuss a possible 
evolutionary scenario for ABDorBa and ABDorBb in terms 
of an earlier quadruple system that evolved to its present state 
via Kozai cycling and tidal friction (Mazeh & Shaham 1979; 
Fabrycky & Tremaine 2007). 

2. Observations and data reduction 

We have carried out three epochs of observations with the 
Long Baseline Array (LBA), the Australian VLBI Network, be¬ 
tween 2007 and 2013 (see Table [T}. Each observation lasted 
12 hours at the frequency of 8.4GHz. Both RCP and LCP po¬ 
larizations were recorded with a rate of 1024 Mbps (2 polar¬ 
izations, 8 subbands per polarization, 16 MHz per subband, 2 
bits per sample), except at Hobart and Ceduna, with a record¬ 
ing rate of 512 Mbps (2 polarizations, 4 subbands per polariza¬ 
tion, 8 MHz per subband, 2 bits per sample). We used the phase- 
reference technique, interleaving scans of the ICRlQ-dcfining 
source BLLacPKS 0516-621 and the star ABDorB (separated 
by 3.6°). The observation sequence target-calibrator-target lasted 
about four minutes. 

We reduced the data using the Astronomical Image 
Processing System (AIPS) program, of the National Radio 
Astronomy Observatory (NRAO), following standard proce¬ 
dures: (i) we calibrated the visibility amplitude using system 
temperatures and antenna gains; (ii) we removed the ionospheric 
contribution (using GPS-based Global Ionospheric Mapo; and 
corrected the parallactic angle; (iii) we performed a fringe- 

1 International Celestial Reference Frame. 
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Table 2. Circular Gaussian fits corresponding to the VLBI maps 
of the components of AB Dor B 


Epoch 

Component 

Flux 

(mJy) 

Diameter 

(mas) 

2007.863 

Ba 

0.82 ±0.18 

2.81 ± 0.10 


Bb 

0.88 ± 0.17 

3.04 ± 0.06 

2010.816 

Ba 

1.39 ±0.08 

2.92 ± 0.02 


Bb 

0.60 ±0.11 

2.02 ±0.13 

2013.625 

Ba 

0.92 ± 0.07 

2.73 ± 0.07 


Bb 

0.63 ± 0.07 

1.74 ±0.02 


search on the calibrator to remove residual contributions to the 
phases; and (iv) we interpolated these solutions from the calibra¬ 
tor onto the star data. 

We imaged the radio sources with the program DIFMAP 
(Shepherd et al. 1994) and cross-checked the results with the 
IMAGR task in AIPS. To re-scale the visibility amplitudes of 
the the calibrator PKS 0516-621, known to be a variable source 
(Sadler et al. 2006; Murphy et al. 2010), we used the 8.4GHz 
flux density values obtained from ATCA measurements of this 
calibrator taken at the same time as our LBA observations (0.85, 
1.40, and 1.50 Jy for epochs 2007.863, 2010.816 and 2013.625, 
respectively). We iterated amplitude and phase self-calibrations 
with deconvolutions using the clean algorithm to finally obtain 
the uniformly-weighted maps of PKS 0516-621 (see Fig.Q]). For 
each epoch, this iterative procedure allowed us to determine both 
the amplitude scaling corrections and self-calibrated phase for 
each telescope. Back to AIPS program, these corrections were 
then interpolated and applied to the AB DorB data. The phase- 
referenced naturally-weighted images of AB Dor B are shown 
in Fig. [2] We emphasize that in this phase-reference mapping 
process the positional information of AB Dor B with respect to 
the external quasar is conserved, thus relating the position of 
AB Dor B to the ICRF. 

3. Results 

3.1. Maps of AB DorBa/Bb 

At each of the three observing epochs, two sources are clearly 
distinguishable in Fig. [2j we identify both spots with the com¬ 
ponents Ba and Bb previously seen in NIR images (Janson et al. 
2007). We determined the flux and position of both components 
from a circular Gaussian least-squares fit (see Table |2}. Both 
components appear unresolved for all the three epochs.These are 
the first VLBI images of AB DorB, where we confirm that both 
components are compact and strong radio emitters. Previous 
studies (Lim 1993; Wolter et al. 2014) already found the system 
to have an intense radio emission with rapid rotation rate (0.38 
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Fig. 1. Contour images of the calibrator PKS 0516-621 at the three LBA epochs. In each map, the lowest contour level corresponds 
to 3 times the statistical root-mean-square (3, 1.5, and 3 mJy beam -1 ) with a scale factor between contiguous contours of V3. The 
peak flux densities in the images are, respectively, 0.89, 1.54, and 1.44 Jy beam -1 . Notice that the second epoch has only intra- 
Australian baselines and therefore a lower resolution. See Table 1. 


days) and strong magnetic activity, as shown by the high X-ray 
luminosity, close to saturation limit for active late-type stars. Our 
radio maps show that each component Ba and Bb should retain 
a high rotation rate to maintain such a intense activity. Actually, 
the AB Dor B IR spectral analysis reported in Wolter et al. (2014) 
is compatible with two rapid rotators (although a model based on 
a single high-rotation star is not excluded). Using the fluxes and 
sizes of Table[2] we derived a mean radio luminosity for Ba and 
Bb of L r = 2.7 x 10 14 erg Hz -1 s _1 , similar to other radio stars 
members of the AB Dor-MG, indicating that gyro-sinchrotron 
emission originated at the stellar corona is the responsible mech¬ 
anism for the radio emission. 

For a proper astrometric analysis, we still need to identify the 
components Ba and Bb, as labeled by Janson et al. (2007) in their 
NIR images, in each of our images of AB DorB. These authors 
selected component Ba to be the brightest one; however, the flux 
density of both stars is very variable at radio wavelengths (see 
Table [2} making an identification under the same criteria diffi¬ 
cult. To break this ambiguity, and based on the similar relative 
position of components Ba and Bb in both the NIR and radio im¬ 
ages, we selected the easternmost component of the radio maps 
to be Ba. This choice was confirmed to be correct after the anal¬ 
ysis of the orbital motion (see Sect. 3.2). It must be emphasized 
that the appropriate registration of Ba and Bb through the three 
epochs is guaranteed by the use of the quasar PKS 0516-621 as 
an external astrometric reference. 


3.2. Orbit determination of AB DorBa/Bb 

To estimate the orbital elements of AB Dor Ba/Bb, we used the 
positions of AB DorBa/Bb resulting from the astrometric anal¬ 
ysis of our three LBA observations, both the relative positions 
(as measured directly on the maps shown in Fig. [2]) and the 
absolute ones (in turn referenced to the position of the quasar 
PKS 0516-621). We included the NIR relative positions avail¬ 
able in the literature (Wolter et al. 2014, and references therein) 
in our fit. We complemented our data sets with the absolute po¬ 
sitions reported in G06 and re-interpreted here. These authors 
made a first attempt to estimate the orbital elements with LBA 
epochs ranging from 1992 to 1996, but their least-squares analy¬ 
sis did not converge to any plausible solution. A possible reason 
for this non-convergence could be the misidentification of the - 


only- component detected, wrongly associated with AB DorBa 
at all epochs. Preliminary fits of the G06 positions along with 
the new ones presented in this paper show that only the posi¬ 
tions at epochs 1992.685 (corresponding to Bb) and 1993.123 
(corresponding to Ba) are compatible with the new LBA data. 
Accordingly, only those positions are included in our fit. The 
rest of the positions of AB Dor Ba/Bb in G06 are therefore not 
included in our analysis. Table [3] summarizes all the positions 
used in this work. 

To find the astrometric and Keplerian parameters of the 
AB Dor B system, we selected a non-redundant data set from the 
positions shown in Tabled In practice, we combined, in our fit, 
the absolute positions of AB Dor Bb and all relative positions of 
AB DorBb to AB DorBa. The reason for this choice is that the 
absolute orbit of AB DorBb and the relative orbit constructed 
using AB DorBa as reference (Bb-Ba) are identical, except for 
their semimajor axes, which are related by the mass ratio. An 
equivalently valid choice would have been a combination of ab¬ 
solute data of Ba and relative data constructed as Ba-Bb. We no¬ 
tice that other combinations of absolute and relative data would 
have produced an extra difference of n radians in the longitude 
of the periastron (oj) between the absolute and relative orbits, 
which would have complicated our orbit analysis (see below). 
We used a least-square fit approach similar to that described in 
G06, but slightly improved to deal with relative data, that we 
outline here. 

We define the absolute position of AB DorBb (aBb,^Bb) at 
epoch t through the expressions: 

g-bAO = a(t 0 ) + li a (t - t 0 ) + Q a (t - f 0 ) 2 + rrP a + 
S a (t.A.B,F,G.P,e,T 0 ) 

<5Bb(0 = b(fo) + Ps(t - to) + Qg(t - to)- + nPg + 

Ss(t,A, B, F, G, P, e, Tq), 

where to is the reference epoch, fi a and /jg are the proper motions 
in each coordinate, Q a and Qg are the secular perspective ac¬ 
celerations in each coordinate (i.e., time variation of the proper 
motions), n is the parallax, P a and Pg are the parallax factors 
(Green 1985), and S a and S „ are the absolute orbital motions of 
AB DorBb in a and 6, respectively. The acceleration terms are 
intended to model the long-term expected curvature of the sky 
trajectory of AB DorB due to the gravitational pull of AB Dor A 
9" away (see G06). We have used the Thieles-Innes elements 
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Fig. 2. Clean maps of the binary ABDorB at the three LBA 
epochs. In each map, the lowest contour level corresponds 
to 3 times the statistical root-mean-square (0.08, 0.19, and 
0.12 mjy beam -1 j with a scale factor between contiguous con¬ 
tours of yfl. The peak flux densities in the images are, respec¬ 
tively, 0.86, 1.25, and 0.89 mjy beam -1 . For image parameters 
see Table [2] In all the maps we have centered the position of 
AB DorBa at the origin. Absolute positions of both Ba and Bb 
are in Table [3] 


(A, B, F, G ) that are combinations of the semimajor axis of the 
absolute orbit a Hh , the inclination i, the longitude of the ascend¬ 
ing node Q, and the longitude of periastron <x>. The three remain¬ 
ing non-linear parameters are the period P, the eccentricity e, 
and the time of periastron passage Tq. 

Similarly, the relative position (ar re i, Ad, constructed as 
Bb-Ba) at epoch t' (not necessarily different from t), are in¬ 
cluded in the fit through the following expressions: 

Aa rc i(T) = Saif, q,A , B, F, G, P, e, T 0 ) 

A6 K [(t') = S s (t', q,A, B, F, G, P, e, T 0 ), 

where q is the ratio between the semimajor axes of the relative 
and absolute orbit, a re i/aBb- We notice that q behaves as a scale 
factor between the absolute and relative orbit, both sharing all 
the orbital elements but the semimajor axes (provided that the 


absolute and relative data sets have been chosen appropriately, as 
explained above). Accordingly, the true Thieles-Innes constants 
of the relative orbit can be written as ( qA , qB, qC , qG), so that 
q remains as the only additional free parameter in our fit, when 
adding the relative orbit data. Given the definition of both data 
types, the x 2 to be minimized is the following: 


X 


z 

i= 1 


(<*Bb0/) ~ &Bb(A)) 2 

°iL iu) 


M 

2 

/= 1 


(A« rc |(t') - Aor re |(f')) 2 


+ Z 

i= 1 
M 

z 

;=i 


O) 

(Ad rd (t;) - Ad re i(t;)) 2 
cr? (f.) 

5rd V i' 


where N is the number of absolute positions, M the number of 
relative positions, the cr’s are the corresponding standard devi¬ 
ations, and the circumflexed quantities are the theoretical po¬ 
sitions derived from the a priori values of the astrometric and 
orbital parameters. We have used the Levenberg-Marquardt al¬ 
gorithm (L-M; e.g.. Press et al. 1992) to minimize the x 2 . As 
other methods to find minima in a non-linear^ 2 space, the effi¬ 
ciency of the L-M algorithm is much improved if good a priori 
values are used for the orbital elements. In our particular case, 
exploration of the relative orbit data yields a reasonable initial 
value for the period: as seen in Table[3j the two NIR positions are 
similar, indicating that the time difference between both epochs, 
nearly one year, could be a promising initial value for the pe¬ 
riod. In addition, the separation and the position angle in the 
two first LBA epochs (2007.863 and 2010.816) are almost coin¬ 
cident, with a difference of less than two milliarcseconds. This 
strongly suggests that the time difference between the two LBA 
epochs (2.95 years) should be approximately an integer num¬ 
ber of the true orbital period. If we assume that three complete 
orbits have elapsed between the two epochs, the orbital period 
would be 0.98 yr, coincident with the estimate made from the 
relative infrared data. Actually, the use of an a priori value for 
the period of 1 yr did facilitate the convergence of the L-M al¬ 
gorithm. Of course the previous argumentation is also valid for 
shorter periods as long as P ~ 1 jn yr, with n being an integer 
number; however, no convergence is found in our fit for periods 
with n > 2. 

The set of orbital parameters that produces a minimum in the 
X 2 (1.2) is shown in Table [4] Plots of the relative and absolute 
orbits are shown in Figs. [3] and [4] respectively. These parame¬ 
ters coincide with those of Wolter et al. (2014) to within their 
standard deviations (except for Q, with a difference of ~ zr/2, 
which may reflect a quadrant ambiguity). The sum of the masses 
of the system can be calculated from a re i and the period P us¬ 
ing the Kepler’s Third Law (/»B a + wbb = a^/P 2 ), which re¬ 
sulted to be 0.53+0.05 M 0 , somewhat lower, but within uncer¬ 
tainties, than Wolter et al. (2014) estimate of 0.69^4 M 0 . More 
importantly, since our data also provide the semimajor axis of 
the absolute orbit of AB DorBb, «b b- we can also calculate the 
mass of the component Ba, ;»B a , using the Kepler’s third law in 
the form m 2 /(rnsa + f«Bb) 2 = a B\,/P 2 - This yields a value of 
m b a =0.28+0.05 M 0 ., with the uncertainty calculed from propa¬ 
gation of the standard deviations of the semi major axis «Bb. the 
period P, and + m^b- Since we have also absolute data of 
AB DorBa, we can estimate, independently of the previous fit, 
a value of m^b (he., not just the difference between mass sum 
and ;«B a )- To do this, we repeated the fit described throughout 
this section but using the absolute positions of AB DorBa and 
the relative positions constructed as Ba-Bb. As explained above, 
this is an equivalent fit to the previous one, except that it provides 
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Table 3. Compilation of astrometric measurements for the AB Dor B system 


Epoch 

Instrument 

Relative positions AB DorBb 
Aa (mas) 

- AB DorBa 

A 6 (mas) 

Reference 

2004.091 

VLT (IR) 

-56.8±3.0 

-33.1±3.0 

(1) 

2005.019 

” 

-64.6±3.0 

-27.4±3.0 

(1) 

2005.909 

” 

-66.7±3.0 

-4.0±3.0 

(1) 

2008.650“ 

” 

9.6±3.0 

-16.4±3.0 

(1) 

2008.855 

” 

-61.3±3.0 

-9.9±3.0 

(1) 

2008.967 

” 

-61.5±3.0 

-24.8±3.0 

(1) 

2009.003 

” 

-57.3±3.0 

-26.7±3.0 

(1) 

2009.131 

” 

-45.6±3.0 

-32.7±3.0 

(1) 

2007.863 

LBA (radio) 

-62.0±0.3 

-10.5±0.7 

(2) 

2010.816 

” 

-60.3±1.5 

-9.7±1.4 

(2) 

2013.625 

” 

-31.3±0.4 

9.3±0.9 

(2) 



Absolute positions AB DorB (LBA) 


Epoch 

Component 

RA (h m s) 

Dec (°'") 


1992.685 

Bb 

5 28 44.41973 ± 0.00060 

-65 26 47.0047 ± 0.0021 

(3) 

1993.123 

Ba 

5 28 44.40441 ± 0.00080 

-65 26 46.9869 ± 0.0028 

(3) 

2007.863 

Ba 

5 28 44.57761 ± 0.00008 

-65 26 45.1002 ± 0.0010 

(2) 


Bb 

5 28 44.56766 ± 0.00008 

-65 26 45.1107 ± 0.0010 

(2) 

2010.816 

Ba 

5 28 44.61098 ± 0.00019 

-65 26 44.7132 ± 0.0008 

(2) 


Bb 

5 28 44.60130 ± 0.00014 

-65 26 44.7229 ± 0.0008 

(2) 

2013.625 

Ba 

5 28 44.63954 ± 0.00015 

-65 26 44.2920 ± 0.0009 

(2) 


Bb 

5 28 44.63453 ± 0.00013 

-65 26 44.2827 ± 0.0008 

(2) 


Notes. a In Wolter et al. (2014) this position was rather considered as an upper bound of 19 mas for the separation of Ba/Bb at 2008.855. (1) Wolter 
et al. (2007); (2) This study; (3) Guirado et al. 2006. The standard deviation of the relative position corresponds to the SNR-based uncertainty 
of the peaks of brightness of ABDorBa and Bb. The absolute position were obtained with reference to the IERS coordinate of the external 
quasar PKS 0516-621 (a = 5 h 16 m 44 s .926178 ,6 = -62° 7' 5".38930). The standard deviation of the absolute position includes, in addition to their 
uncertainty of their respective peak of brightness, the contribution of the propagation media and the reference source structure. 


an estimate of a Ba instead of «Bb- Again, using the Kepler’s Third 
Law, now with a Bd , we obtain a value of mBb=0.25+0.05 M 0 , be¬ 
ing this estimate largely independent of m Bd and the total mass. 

Regarding the astrometric parameters in Table |4j we notice 
that the long-term orbit of the pair AB Dor B around AB Dor A 
is reflected on the values of the perspective acceleration; actu¬ 
ally, solving for Q a and Q n reduced the rms of the residuals 
by a factor of three. The magnitude of the acceleration can be 
easily reproduced from the simple expression (circular orbit) 
(2n/P AB ) 2 x tab. taking the A/B distance r AB = 9.22" at the 
reference epoch (2000.0) and assuming reasonable values of the 
period of B around A, P AB (~2000 yr; within the plausible orbits 
defined in G06). 

Finally, as seen in Figs. [3] and [4j the orbit of the pair in 
AB DorB is not fully covered. However, the orbital elements in 
Table [4] are satisfactorily determined. This apparent contradic¬ 
tion may be explained from the following two facts. First, the 
so-called “Eggen’s effect” (Lucy 2014), which predicts that for 
poorly covered orbits, the quantity a 2 /P 2 has little variation for 
values of the parameters (a, P) that minimize the x 2 defined in 
Sect. 3.2 (note that a 2, /P 2 is proportional to the star masses, and 
therefore sensitive to any trace of orbital motion evidenced in 
one of the components of the binary). And second, the fact that 
we have a precise a priori value for the orbital period, which 
combined with the previous effect, leads to a firm estimate of 
the semimajor axis, greatly constraining the rest of the orbital 
elements. 

4. Discussion 

4. 1. Stellar evolution models for PMS stars 

The measurements of dynamical masses are essential to check 
PMS stellar evolution models. The mass is the most important 


stellar parameter, but determinations of the luminosity and the 
temperature are also necessary to calibrate the theoretical mod¬ 
els. Extensive description and comparison between existing stel¬ 
lar models can be seen in Hillenbrand & White (2004). 

In this paper, we have considered isochrones and isomasses 
corresponding to the PMS models of Baraffe et al. (1998; 
BCAH98), Siess et al. (2000; S00), Montalban & D’Antona 
(2006; MD06), and Tognelli et al. (2011, 2012; TDP12). A 
metallicity value of [Fe/H |= 0.0 has been adopted as the av¬ 
erage value of the AB Doradus moving group has reported to be 
[Fe/H |=0.02 + 0.02 (Barenfeld et al. 2013). All the quoted mod¬ 
els provide theoretical masses for stars down to 0.1 M 0 , but for 
different mass spacing. Low-convection efficiency (i.e. mixing 
length parameter otml = 1 -0), which optimizes the comparison 
with our data, have been used in the case of BCAH98, MD06, 
and TDP12 models. For S00 models, tracks are available only 
for solar calibration, with uml = 1.61. 

Besides the convection treatment, other differences among 
the quoted PMS low mass star models are relevant, in particu¬ 
lar, the adopted equation of state, the radiative opacity, the at¬ 
mospheric structures used to specify the outer boundary con¬ 
ditions, and the adopted initial abundances of the chemical el¬ 
ements (mainly helium and total metallicity) at a fixed [Fe/H] 
value (see, e.g., Siess et al. 2001, Mathieu et al. 2007, Tognelli 
et al. 2011). As a result, the age and mass inferred by compar¬ 
ing evolutionary tracks and data might significantly change for 
different stellar models (i.e. input physics). 

These models are displayed in Fig. 0 where different iso¬ 
masses and isochrones are plotted. We have placed AB DorBa 
and Bb in the HR diagrams of Fig. fusing the spectroscopically 
determined effective temperatures of Wolter et al. (2014) and the 
K magnitude measured by Janson et al. (2007). BCAH98 and 
TDP12 models suggest that AB DorBa and Bb are coeval, with 
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Fig. 3. Relative orbit for the binary ABDorBb using to the orbital elements in Table 4 (with a le i). The ABDorBa component is 
indicated by the asterisk at the origin. Star symbols and epochs correspond to the VLBI data. For the sake of clarity, the NIR points 
are marked with numbers, following a chronological order which corresponds with their entries in Table 3. A prediction of the 
relative position at epoch 2008.605 is shown as an empty circle, indicating that the upper bound of 19 mas suggested by Wolter et 
al. (2014) is fully accomplished. 



Fig. 4. Absolute orbits of ABDorBa and ABDorBb using the orbital elements in Table 4 (with « H a and asb- respectively). The 
positions of components Ba (triangles) and Bb (star symbols) are marked along with their respective epochs. The center of mass 
(CM) of the system is placed at the origin. 


an age between ~50 Myr and lOOMyr, similar to that reported 
by Janson et al. (2007) and Wolter et al. (2004). S00 provides 
a slightly older age, but non-coevality. Finally, MD06 predicts 
substantially older ages, 100-125 Myr for both stars. The latter 
age range seems to agree with the age derived by Barenfeld et 
al. (2013) for the nucleus of the AB Dor-MG (>1 lOMyr). 


On the other hand, the theoretical masses predicted from the 
models agree with our dynamical estimates just at the extreme 
of their uncertainties. The models of BCAH98 and TDP12 (with 
a denser distribution of isomasses than S00) predict mass ranges 
of 0.18-0.25 M 0 for component Ba and 0.14-0.21 M 0 for com¬ 
ponent Bb. MD06 models also underpredict the masses of Ba 


6 









Azulay et al.: Dynamical masses of the binary low-mass star AB Doradus B 



Teff(K) 


- 

1 

0.30 M 0 

logft) — 7.0 

0.20 M 0 


°- 10 M o ’ 

■ 

- 

■ 


log(t) = 8.0 

^4 

Siess et al. 2000 



3600 3400 3200 3000 


Teff(K) 




Fig. 5. Comparison of ABDorB components with some PMS theoretical models (Baraffe et al. 1998, top left ; Siess et al. 2000, 
top right ; Montalan & D’Antona 2006, bottom left; Tognelli et al. 2011; bottom right). For each model isomasses (solid lines) and 
isochrones (dashed lines) are plotted. We highlight the nearest tracks available corresponding to our dynamical mass values. The 
theoretical masses are consistent with our dynamical estimates just at the extreme of their uncertainties. 


and Bb but with a better agreement, in practice within the stan¬ 
dard deviations of our dynamical measurements. Therefore, for 
the BCAH98, TDP12, and S00 models, the dynamical masses 
reported in this paper are 30% and 40% larger than the cen¬ 
tral values of the range predicted for AB Dor Ba and Bb, respec¬ 
tively. This percentage of disagreement decreases to 10% for the 
MD06 models. Our results seem consistent with similar com¬ 
parisons of dynamical and theoretical masses done by other au¬ 
thors (Hillenbrand & White 2004; Stassun et al. 2004; Mathieu 
et al. 2007), who concluded that models underpredict stellar 
masses by 10-30% for PMS stars with masses in the range 0.3- 
1.2 M g . Our comparisons suggest that this disagreement holds 
for masses below 0.3 M 0 . We notice that other stellar models de¬ 
signed specifically for lower mass objects (e.g. DUSTY models, 
Chabrier et al. 2000) do not cover the range 0.2-0.3 M e . 


4.2. An alternative scenario for the binary AB DorB 

As noted in Sect. 3.1, both AB DorBa and AB DorBb are strong, 
compact radio emitters, which contributes to the emerging per¬ 
ception that many radio stars reside in double or multiple sys¬ 
tems (Melis et al 2012; Azulay et al. 2014). Considering the 
radio/X-ray activity correlation (Giidel et al. 1993), this percep¬ 
tion finds support in other, better reported statistics which relates 
X-ray activity and binarity, namely, that about 79% of the X- 
ray emitters are in binary systems (Makarov & Eggleton 2009). 
Many of these binaries correspond to close binaries, but interest¬ 
ingly, a good fraction of the total X-ray binaries (67%; Makarov 
& Eggleton 2009) are part of hierarchical triple systems, with 
the third body in a wider orbit. In this scheme, the presence 
of the wide companion induces dynamical perturbations to the 
orbit of the close binary at each periastron passage, through 
the Kozai cycles and tidal friction mechanism (e.g., Mazeh & 
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Table 4. Estimates of the astrometric and orbital parameters of 
AB DorB ab 


Parameter Value 


ffo (h m s) 

5 28 44.48396 ± 0.00022 

So (° ’ ”) 

-65 26 46.0573 ± 0.0013 

Pc (S yr 1 ) 

0.01054 ± 0.00012 

ps (arcsec yr 1 ) 

0.1287 ± 0.0005 

Qa (s yr 2 ) 

0.000008 ± 0.000001 

Qs (arcsec yr 2 ) 

-0.00010 ± 0.00005 

7i (arcsec) c 

0.0664 ± 0.0005 

P (yr) 

0.986 ± 0.008 

Ufcl ( ) 

0.052 ± 0.002 

a ( ) 

0.028 ± 0.002 

°Bb (”) 

0.025 ± 0.002 

e 

0.6 ±0.1 

i(°) 

121 ± 5 

w B b (°) d 

54 ± 20 

DO 

270 ± 15 

To 

2003.68 ± 0.05 

m a (Mo) 

0.28 ± 0.05 

m Bb (Mo) 

0.25 ± 0.05 


Notes. a The reference epoch is 2000.0. b The number of degrees of 
freedom of the fit is 14; the minimum value found for the reduced x 2 
is 1.2. c We note that our parallax estimate is more accurate than the 
Hippcircos value given for AB Dor A and still compatible. d For the ab¬ 
solute orbit of AB Dor Ba m Ba = w Bb + tt. 


Shaham 1979; Fabrycky & Tremaine 2007). Provided that the 
inner and outer orbits are misaligned, such an interaction pro¬ 
gressively produces a loss of angular momentum that leads the 
close binary to shrink, to increment its orbital period, and even¬ 
tually to a merging process. The resulting merger will retain part 
of the angular momentum showing, as a consequence, an ex¬ 
traordinary high rotational velocity, which will translate to high 
levels of activity both at X-ray and radio. This appears precisely 
to be the scenario found in some of the young radio stars in bi¬ 
nary systems. In fact, Makarov & Eggleton (2009) explained the 
origin of the high rotation rate of AB Dor A (0.5 days) via Kozai 
cycling and tidal interaction between a close binary (now merged 
in AB Dor A) and the 0.090 M 0 low mass companion AB DorC, 
orbiting AB Dor A at a mean distance of 2 A.U. with a period of 
11.74 yr. 

Could a similar mechanism be acting in AB DorB? The re¬ 
markable high radio emission in both AB Dor Ba and AB Dor Bb 
reported in this paper and the rapid rotation suspected to occur 
in the two stars, <0.88 days (Wolter et al. 2014), provide sup¬ 
port to this hypothesis. In addition, the separation Ba-Bb is only 
0.15 A.U. at periastron with a much faster period than the pair 
AB Dor A/C, ~ 1 yr. Therefore, a more efficient Kozai cycling 
and tidal interaction might be taking place. Interestingly, this 
hypothesis would imply a double-double scenario where both 
components, Ba or Bb, resulting from their respective mergers, 
would act as mutual "third bodies”, one with respect to each 
other. This double-double scenario appears not to be rare and 
it has been proposed to explain the light curves of quadruple 
systems of eclipsing binaries (Cagas & Pejcha 2012). 

5. Conclusions 

We have shown the first VLBI images of the binary AB Dor B, 
where we detect the presence of compact radio emission in both 


components. The scientific output of ourLBA astrometric mon¬ 
itoring has been optimized, since both the absolute and relative 
orbits have been determined in combination with published NIR 
relative positions. The dynamical masses of the individual com¬ 
ponents of AB Dor B are very similar (0.28 and 0.25 M Q , for Ba 
and Bb, respectively). These values are 10-30% (10-40%), de¬ 
pending of models, larger than the theoretical estimates of PMS 
evolutionary models for AB DorBa (AB DorBb), underlying the 
known tendency of these models to underpredict the masses, 
yet within 2-cr of the predicted values. Comparisons in H-R di¬ 
agrams favour an age between 50Myr and lOOMyr, although 
Montalban & D’Antona (2006) models predicts ages older than 
lOOMyr. ABDorB is one of the systems in the ABDor-MG 
whose two components are radio emitters, contributing to the ev¬ 
idence that many young, wide binaries are strong radio emitters 
(as AB Dor A and HD 160934; see Azulay et al. 2014). This ra- 
dio/binarity correlation may be the radio counterpart of the better 
established relationship between wide binaries and X-ray activ¬ 
ity, which can be explained under the assumption that the present 
wide binaries were originated by dynamical interaction of hier¬ 
archical triple systems. According to this, ABDorBa/Bb could 
be understood as two earlier close binaries, a double-double sys¬ 
tem, where both binaries were forced to shrink and merge via 
mutually-induced, cyclic dynamical perturbations. The evolved 
two mergers (the present stars Ba and Bb) would thus keep part 
of the angular momentum as a very rapid rotation, responsible 
in turn of both the X-ray and radio activity. 

New, more precise estimates of the dynamical masses would 
be desirable via an improved coverage of the orbit, where the res¬ 
olution and sensitivity of the LBA array has shown to be essen¬ 
tial. Likewise, as pointed out by Janson et al. (2007) and Wolter 
et al. (2014), high-resolution spectroscopy of the components in 
AB Dor B would also constrain better their placement in H-R di¬ 
agrams. The latter observations would also help to confirm the 
suspected rapid rotation of both components, as well as to de¬ 
termine the tilt of the spin axis of each star with respect to the 
orbital plane, which could add arguments in favor of the interpre¬ 
tation of the binary AB DorB as an earlier double close-binary 
system. 
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